Coupling of a Single Diamond Nanocrystal to a Whispering-Gallery Microcavity: 
Photon Transportation Benefitting from Rayleigh Scattering 
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We study the Rayleigh scattering induced by a diamond nanocrystal in a whispering-gallery- 
microcavity-waveguide coupling system, and find that it plays a significant role in the photon 
transportation. On one hand, this study provides a new insight into future solid-state cavity quan- 
tum electrodynamics toward strong coupling physics. On the other hand, benefitting from this 
Rayleigh scattering, novel photon transportation such as dipole induced transparency and strong 
photon antibunching can occur simultaneously. As potential applications, this system can function 
as high-efficiency photon turnstiles, fn contrast to [B. Dayan et at, Science 319,f062 (2008)], the 
photon turnstiles proposed here are highly immune to nanocrystal's azimuthal position. 

PACS numbers: 78.35.+C, 42.50.Pq, 42.50.Ar 



Cavity quantum electrodynamics (CQED) studies 
light-matter interactions inside a resonator, which offers 
an ideal platform for quantum optics 1]. Tremendous 
progress has been made by coupling single dipoles to 
different microcavities (for a review, see @). Among 
them, whispering-gallery mode (WGM)-type microcav- 
ity [1] is most promising due to its extremely high Q fac- 
tor, small mode volume, excellent scalability and ease for 
low-loss transport of nonclassical states using an optical 
fiber [3, IK . On the other hand, nitrogen- vacancy (NV) 
centers have recently emerged as an important candidate 
for quantum information processing because they pos- 
sess long-lived spin triplets at room temperature § 0] ■ 
Combining high-Q WGM microcavities and NV centers 
represents a promising solid-state CQED system, and at- 
tracts much attention very recently [4, 5;, Q ■ 

One of the distinct properties of WGMs is that they 
belong to traveling waves, unlike standing modes in a 
Fabre-Perot cavity. In other words, WGM microcavi- 
ties typically support twin modes, clockwise (cw) and 
counter-clockwise (ccw) propagating waves with a degen- 
erate frequency. In this paper, we investigate the inter- 
action of the twin WGMs coupled to a NV center in a 
diamond nanocrystal. On one hand, we find that not 
only the dipole transition of the NV center but also the 
Rayleigh scattering by the nanocrystal itself play signif- 
icant roles in the coupled system. Thus strong coupling 
condition in such CQED system should be re-defined. On 
the other hand, nonclassical effects are predicted bene- 
fitting from the nanocrystal scattering. Dipole induced 
transparency (DIT) Q appears in the presence of the 
strong Rayleigh scattering, accompanying strong pho- 
ton antibunching in transmitted fields. Moreover, high- 
efficiency photon turnstiles can be implemented highly 
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FIG. I: Schematic illustration of the nanocrystal-microcavity- 
waveguide coupling system. 



immune to the azimuthal position of the nanocrystal. 

Figure [T] describes the present coupling system. A di- 
amond nanocrystal embodying a single NV center is lo- 
cated on the surface of a microtoroidal cavity [l(J. A 
tapered fiber is used to couple light with frequency lo p 
in and out of the microcavity. Under the rotating wave 
approximation, the Hamiltonian of the coupled system 
can be written as (H = 1) 

H = Hq + H i + H2 , 
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Here summation indices (m, m!) run through cw, ccw; a m 
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and bkj (k — 1,2,3) denote annihilation operators of the 
cavity and the reservoir modes, respectively; <r_ = \g) (e| 
and a+ = |e) {g\ stand for the descending and ascending 
operators of the dipole transition at zero-phonon line. 
The Hamiltonian Ho describes the free evolution of the 
system consisting of the NV center, WGMs and reservoir 
modes, with frequency w e , lo c , uj^j, respectively. The 
first (second) term in H\ characterizes the dipole in- 
teractions between the NV center and WGMs (reservoir 
modes), with the coupling strength G m (Gj). The first 
term in H 2 describes the nanocrystal-induced scatter- 
ing into the same (to = m') or the counter-propagating 
(to 7^ to') WGM fields with strengths g m ,m' > while the 
second term represents the WGMs-reservoir scattering 
coefficients g m ,j- 



In the full quantum theory, the coherent cou- 
pling strengths G m and Gj can be calculated by 

G m = G = ^{u Jc /{2he () e s V c )) 1/2 f c {^) and Gj = 

li{uj lj /{2he e s V 1 j)) 1/2 , where /z = 2.74 x 1CT 29 C-m 
represents the dipole moment of the NV center transi- 
tion. So is the electric permittivity of the vacuum and 
e s = 1 denotes the relative permittivity of the surround- 
ing medium (vacuum here). V c and V\j are the quantized 
volumes of the cavity and the reservoir modes interact- 
ing with the dipole, respectively. f c (f) = \E(r)/E max \ 
is the normalized field distribution function of WGMs. 
For the subwavelength nanocrystal, the above scatter- 
ing interaction can be modeled in a dipole approxima- 
tion where the electric field of the input wave (either 
WGMs or reservoir modes) induces a dipole moment in 
the scatterer, namely, scattering induced dipole. In the 
case of elastic Rayleigh scattering, the coupling strengths 
are g m>m > = -g = -af%(r)u c /(2V c ) and g m ,j = ~9r = 
-af c {f)uj c (n m • n/ CJ )/(4V c V fcj ) 1 / 2 (m = cw, k = 2 or 
m = cw, k = 2)[lll- Here a = 47ri? 3 (e d - e s )/(e d + 2e s ) 
is the polarizability for the spherical nanocrystal with 
radius R where Ed — 2.4 2 denotes the electric permittiv- 
ity of the diamond. fi m and fik.j are unit vectors of the 
fields. 



We can find that both the NV center-reservoir and 
WGM-reservoir coupling interactions have the same for- 
mat. By using Weisskopf-Wigner approximation, the 
coupling of the NV center to the reservoir can be re- 
garded as a decay of the excited state, i.e., the well- 
known spontaneous emission with rate "f s ; while the scat- 
tering of WGMs to the reservoir can be modeled by 
an energy damping of WGMs, with the damping rate 

kr = a 2 f^(f)s^ 2 U!^/(6Trc 3 V c ) [ll|, where c is the speed 
of light in vacuum. Transforming travelling to standing 
modes, a± — (a cw ± a ccw )/V2, the equations of motion 



for the coupled system are given by [12J 

^± = (2*. 9 - K+ )a + - sftiGo. - MaZ + f+, (2) 



da,- 
~dt 
do~- 
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(3) 



jU =-{iA ec +^-)a-+iV2Ga+a z +$ 1 . (4) 

Here a z = |e) (e\ ~ \g) (g\, A oc = u) e - oj c , k + = 
kr + {ko + K i)/2, K- — {ko + Ki)/2; ko = u c /Q denotes 
the intrinsic damping of the WGMs with Qo being the 
intrinsic quality factor; k\ stands for the cavity-taper 
coupling strength; a 1 ^ is the input field. The opera- 
tors f + , f_, fi are the noise operators that conserve the 
commutation relations at all times. One can find that 
the symmetric mode a+ is strongly coupled to the NV 
center (through dipole interaction) and the nanocrystal 
(through the scattering), while the anti-symmetric mode 
a_ keeps uncoupled to them. 

When the cavity is excited by a weak monochromatic 
field (i.e., weak- field approximation), the NV center is 
predominantly in the ground state. Thus cr z can be sub- 
stituted for its average value of —1, and Eq. ([4| becomes 
linear. Utilizing the standard input-output formalism 
Kio ra [l3|, we obtain the outputs 
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where D + = i(A pc + 2g) - n + + 2G 2 /[i(A pc - A cc ) - 
7 s /2], D_ = zA pc - K-, with A pc = u p — w c , % 
and fg are noise operators. We take the cold reser- 
voir limit where the reservoir modes are all initially 
in vacuum states. Thus the expectation values of all 
noise operators can be neglected because they are an- 
nihilated when acting on the initial vacuum states. Then 
the transmission (a°^a°^) / (aj^aj. 1 ^) and reflection 

(a™4 ta ccw) / (4wacw) can be determined. 

Under practical parameters j s /2ir = 13 MHz, Qo 
= 10 8 , m = 0.5k , V c - 200 ^m 3 , f c (?) - 0.47, we 
obtain {G, n , /2ir - {180, 4.7, 2.35}MHz. We will 
use these parameters unless otherwise specified. Figure 
[2fa) plots these parameters depending on the particle ra- 
dius R. We can approximately divide Fig. into three 
regions. In region I, G far exceeds any other parameters 
for a small diamond nanocrystal. The NV center couples 
to the twin WGMs, creating a pair of standing modes. 
There exist three separate resonance dips in transmission 
spectrum: one at the original zero detuning due to the 
decoupling between the anti-symmetric mode and the NV 
center, and the other two splitted by 2\[2G due to the 
strong coupling between the symmetric standing mode 
and the NV center. The linewidth of the central dip is 
Ko + K\ ; while the linewidths of the two sideward dips are 
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FIG. 2: (a) Parameters {G, g, j 3 , kq, Ki, kb} vs. the radius R 
of the nanocrystal. (b) The transmission spectra for various 
R. Here Q = 10 8 , m = 0.5k and A ec = 0. 



kr + {kq + K± + 7s)/2, which can be interpreted by the 
dressed state generated by the strong coupling. 

With the increase of R, the coefficient g grows rapidly. 
In region II, g becomes comparable with G, while kr is 
still much smaller than G because g(nR) oc R 3 (R 6 ). The 
scattering not only moves the two sideward resonances 
with detunings —g ± (g 2 + 2G 2 ) 1 / 2 , but also increases 
their linewidths and decreases the coupling efficiencies. 
When the size of the nanocrystal is large enough (region 
III), the scattering damping rate kr becomes comparable 
to or even exceeds G. In this situation, the scattering 
damping rate is so large that the two sideward dips nearly 
vanish. In all cases, the central dip keeps unchanged, and 
its coupling depends mainly on kq and K\. 

Transmission spectra for different R varying from 10 
to 70 nm are ploted in Fig. [2fb) , confirming the above 
analysis. As the nanocrystal-induced scattering does 
occur, the new strong coupling condition requires that 
G ^> {k + ,j s }, corresponding to regions I and II. This 
provides new insights into future solid-state CQED. 

In the region where the present strong coupling condi- 
tion does not match, interesting DIT effects in the trans- 
mission spectrum can be predicted. Let us consider the 
conditions for a general DIT phenomenon. First, in the 
absence of the dipole, the new critical coupling condi- 



tion (k{ 



Ag ) is required to create a relatively 



broad platform with a low transmission. This can be 
implemented though the destructive interference among 
the directly transmitted field uncoupled to the microcav- 
ity and the output field via cavity modes, as discussed 
in detail in the following. Second, the linewidth of the 
resonance should be larger than the dipole-cavity cou- 
pling strength > G, i.e., the bad cavity limit). 
These conditions are verified in Fig. 3, which presents 
the transmission, reflection and energy loss in different 
conditions. 

In the weak scattering case (g -C Ko), the initial crit- 
ical coupling condition without dipole is simplified as 
Ki = kq. For cavities with high intrinsic quality fac- 
tor, the mode splitting occurs instead of DIT, as shown 
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FIG. 3: Transmission, reflection and energy loss of the present 
system working on the critical coupling point Ki = (kq + 
4g2-ji/2 (a)-(d), weak scattering case with R = 6 nm and 
A cc = 0. (e)-(h), strong scattering case with R = 60 nm and 
A cc = -g - G 2 /g. The high and low Q are 10 8 and 10 6 , 
respectively. 



in Figs. 3(a)-3(b). This is because the dipole-cavity cou- 
pling strength exceeds the resonance linewidth. While a 
cavity with low intrinsic quality factor can produce DIT 
as shown in Figs. 3(c)-3(d), but there are significant en- 
ergy losses, yielding the low DIT peak, as well as low 
on/off contrast ratio. 

In the strong scattering case (g 3> kq), however, we 
find that the nanocrystal-induced Rayleigh scattering 
plays a constructive role in generating DIT and can sig- 
nificantly enhance it, as demonstrated in Figs. 3(e)-3(h). 
In general, the field radiated by the NV center dipole 
takes part in the interference mentioned above. If we 
control the detuning A oc = —g — G 2 / g so that the radi- 
ated field intensity has a maximum value at A pc ~ — g 
where the initial critical-coupling point is located at. The 
DIT window (also at A pc ~ —g) with a high efficiency 
(the peak value exceeds 0.8) appears, because the intrin- 
sic loss is far smaller than the cavity-taper coupling rate 
(«i - 2g > k ). 

High-efficiency PIT holds great potential for pho- 
ton turnstiles |13l . flij . We are thus interested in 
the photon statistics of the output [n| [3, de- 

(2) 

scribed by the second-order correlation functions gin = 
((a™ t t) 2 (a™ t ) 2 )/(a™ tt a™ t ) 2 . Figure Ha) depicts the 
strong photon antibunching of the transmitted field with 
a large bandwidth, corresponding to the case in Fig. Elf)- 
This is because the net transmitted field mainly orig- 
inates from the NV center dipole and a single dipole 
cannot emit two or more photons simultaneously. Ac- 
tually, excitation of the system by a first photon blocks 
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FIG. 4: Normalized second-order autocorrelation functions 
Sew (a) and phase shift (b) of the corresponding transmitted 
field in Fig. 3(f). The dashed vertical lines and blue arrows 
indicate the position of A pc = —g. Inset: Four components of 
the transmitted field in complex plane (see text). 

the transmission of a second photon, known as photon 
blockade [I?} • It is of importance that the strong anti- 

(2) 

bunching (g C w = 0) occurs exactly at the DIT window, 
which guarantees the high efficiency of the antibunched 
light. 

In the DIT window, the phase shift CW of the trans- 
mitted field (defined as 4> cw = a,vg[a°^ / a™ w \) has also 
been plotted in Fig. Hfb), which exhibits a very large 
dispersion in the region near A pc = —g. This enables 
us to control the phase shift and group delay of the out- 
put field. It is interesting that the transmitted field ex- 
periences a 7r/2 phase shift at the DIT peak, different 
from the previous DIT in Q (zero shift). To explicitly 
explain this, from Eq. ([5]), we note that the transmit- 
ted field is an interference of four components: directly 
transmitted field uncoupled to the microcavity, output 
field via the anti-symmetric mode, output field via the 
Rayleigh scattering, and output field via the NV-center 
dipole interaction, denoted by afc=i,2,3,4, respectively. In 
the inset of Fig. 0|b) we plot the four components for 
A pc = — g in the complex plane. Here, we find a\ = a"^ 
and d2 = (Ki/2)aJ.^/(iA pc — «;_), which cannot be tuned 
via the nanocrystal; a 3 = (Ki/2)a c ^,/(z(A pc + 2g) — k + ) 
depends on the scattering strength g while 04 depends 



mainly on G and A cc . It can be seen that a±, ai, a 3 de- 
structively interference, creating zero net field. 04 has a 
large amplitude with the relative phase 7r/2. This further 
explains the strong antibunching of the transmitted field 
in the DIT window. 

Finally, in Ref. [IH, photon turnstiles are realized 
by coupling microcavity to single cooled caesium atoms. 
Standing wave modes built there are induced by scat- 
tering from defects and surface roughness of the micro- 
toroidal cavity, and thus photon turnstiles are sensitive 
to the atom's azimuthal position. Averaging over the 
azimuthal angle results in a significant reduction of the 
on/off contrast ratio of transmission in the DIT window 
(see Fig. 1(E) in Ref. [Hj]). While in the present paper, 
Rayleigh scattering is induced by the diamond nanocrys- 
tal itself. Thus these devices are highly immune to az- 
imuthal position of the nanocrystal. This is of special im- 
portance when multi-nanocrystal are involved in the mi- 
crocavity system. For example, the dynamics of the cou- 
pling system will strongly depend on the inter-positions 
of these nanocrystals. 

In summary, we have analyzed a single diamond 
nanocrystal with a NV center coupled to high-Q counter- 
propagating twin WGMs of a microtoroidal cavity. It 
is found that the Rayleigh scattering induced by the 
nanocrystal itself not only induces strong interaction be- 
tween the twin WGMs, but also results in the scatter- 
ing loss of the WGMs. This provides a new insight into 
future solid-state CQED. We also reveal that Rayleigh 
scattering can play a positive role in obtaining DIT with 
high efficiency, accompanying strong photon antibunch- 
ing. The system can perform as photon turnstiles which 
are highly immune to the azimuthal position of nanocrys- 
tal. Our investigation can be of practical use in controlled 
interactions of single quanta and scalable quantum logic. 

Note added - We note a recent publication [IS] after 
the submission of our work. 
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